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Abstract 



Quantum field theory on the noncommutative two-dimensional Min- 
kowski space with Grosse- Wulkenhaar potential is discussed in two ways: 
Mh, In terms of a continuous set of generalised eigenfunctions of the wave op- 

erator, and directly in position space. In both settings, we find a new 
type of divergence in planar graphs. It is present at and above the self- 
dual point. This new kind of divergence might make the construction of 
Cn ' a Minkowski space version of the Grosse- Wulkenhaar model impossible. 

> 

[^ ■ 1 Introduction 

^^ ■ The extensive study of noncommutative quantum field theories (NCQFT) that 

•Q I started about fifteen years ago has undergone several twists since then (for 

motivations and an overview, we refer to the reviews [T]). Right from the start, 
different approaches were followed. While Filk derived modified Feynman rules 
for the Euclidean case [5], Doplichcr, Fredenhagen and Roberts started on the 
noncommutative Minkowski space and proposed a Hamiltonian setting for the 
S^ ■ quantisation of field theories [3]. In the context of Filk's Feynman rules, the 

i-j I phenomenon of UV/IR mixing [4] was found. For a couple of years, this was an 

obstacle for a systematic treatment of renormalisation. 

The Hamiltonian approach and the modified Feynman rules are equivalent 
when time commutes with all spatial coordinates. However, when this is not 
the case, the naive application of Filk's Feynman rules to field theories on the 
noncommutative Minkowski space leads to a violation of unitarity, in the sense 
that the cutting rules no longer hold [5]. As Bahns et al. |S] pointed out, 
this violation of unitarity is the consequence of an inappropriate definition of 
time-ordering and not present in the Hamiltonian setting [7]. As another al- 
ternative, they proposed the Yang-Feldman approach [5]. Then, the UV/IR 
mixing manifests itself as a distortion of the dispersion relations [9l [10] . 

In recent years, the study of models with an added harmonic oscillator poten- 
tial became popular. Such a modification was first proposed by Langmann and 



Szabo [TT] , who showed that in such a way the action of the complex 0^ theory 
can be made invariant under Fourier transformation. This occurs for a partic- 
ular value of the harmonic oscillator frequency, the so-called self-dual pointj^. 
For the case where the harmonic oscillator potential is replaced by a constant 
magnetic field (of the same frequency), the model is solvable (but trivial) |12| . 
Further evidence for the need for an added harmonic oscillator potential came 
from the success of the Grosse-Wulkenhaar model. They showed that with such 
a modification, the noncommutative (j>'^ model is renormalisable not only in 
two [13) . but also in four spacetime dimensions |14j . Even better, the model is 
asymptotically safe (but not free!), since the /3-function is bounded [ISlfTB] . 

In this approach, one uses (Weyl symbols of) ket-bras of harmonic oscillator 
eigenfunctions to transform the model to matrix form. Then, the interaction 
term takes a particularly simple form. It is precisely at the self-dual point 
that also the propagator becomes simple (diagonal). It turns out that in four 
spacetime dimension^j, the self-dual point is a fixed point of the theory [T21 HB] . 

While the study of Euclidean models with an added harmonic oscillator 
potential was a spectacular success, very little is known about the corresponding 
Minkowski space versionfO. As we show here, the self-dual point is a special 
point also on Minkowski space, but in an unexpected way: It is at this point 
that a strange kind of divergence appears in planar graphs. At first sight, this is 
paradoxical: Usually the planar graphs are as in the commutative case. But in 
two dimensions all scalar field theories are superrenormalisable. Since the degree 
of the singularity of the propagator only depends on the kinetic terrrQ, this is 
also true if a non-constant potential is added. Hence, there is no divergence in 
the commutative case and we would not expect to find one in the planar sector. 
However, the fact that planar graphs are exactly as in the commutative case 
relies on the cancellation of the twisting factors. This, in turn, happens due 
to momentum conservation. But momentum is not conserved because of the 
quadratic potential. Thus, there is no reason to expect that the planar part is 
exactly as in the commutative case. 

We will find these divergences in two different ways. In the first approach, we 
restrict ourselves to the self-dual point and use a continuous set of generalised 
eigenfunctions of the wave operator as a basis for quantisation. As a first step, 
we assume naive Feynman rules and compute the fish graph in the (p^ model. 
While there is some ambiguity stemming from different possible choices for 
the Feynman propagator, one generically finds a peculiar type of divergence in 
planar graphs. It is no UV divergence in that it does not stem from a divergent 
loop integral. Instead, it comes from the kinematical factors at the vertices. If 
the divergence is formally removed, one also finds a violation of unitarity, in the 
sense that the cutting rules are not fulfilled. This is not surprising given that 
we postulated Feynman rules without caring for correct time-ordering. This 



^In this work, this term does only refer to the frequency of the harmonic oscillator potential. 
The models discussed here are in general not self-dual in the sense of |11| . 

^In two spacetime dimensions, the quadratic potential is only needed in an intermediate 
step. In the limit where the cutoff in the matrix base is removed, its frequency vanishes 1131 . 

^In 1171 . first steps in this direction were taken. The relation of their approach to the ones 
discussed here is clarified in Appendix [Al 

*This is the case for the retarded propagator |18l . That the same is also true for the 
Feynman propagator and the two-point function should rather be viewed as a condition for 
suitable states. Basically this is the Hadamard, or in modern terms, the microlocal spectrum 
condition for QFT on curved spacetimes | 19| . 



problem can be cured by quantising a la Yang-Feldman. However, we will argue 
that this does not remove the divergences. 

The second approach that leads to the same conclusion is a treatment of the 
model in position space. In that case, one is not restricted to the self-dual point. 
In order to circumvent the ambiguities connected to a choice of the Feynman 
propagator, we start by explicitly constructing the retarded propagator. It turns 
out to be increasing like a Gaussian in some directions, so we interpret it as a 
distribution on a Gelfand-Shilov space. We show that for distributions on that 
space, the planar ^-product at different points can not be defined via duality 
if one is at or above the self-dual point. Again, this shows that the problem is 
not a UV divergence, since it occurs before taking the limit of coinciding points. 
We also show that when the planar 7k-product is calculated in a formal sense, 
one finds a geometric series that diverges at and above the self-dual point. 

The appearance of this new kind of divergence is an interesting phenomenon, 
that, to our opinion, deserves more detailed studies. In particular, it should be 
checked whether it also occurs in four spacetime dimensions, since there the 
self-dual point is a fixed point of the theory [TSJ [H] , at least in the Euclidean 
case. 

This paper is organised as follows: In the next section, we fix some notation. 
In Section |3] we introduce the continuous basis of generalised eigenfunctions. 
This is used in Section |4] to postulate naive Feynman rules and discuss the fish 
graph of the (f)^ model in this setting. We find the above mentioned violation of 
unitarity and the new type of divergence in the planar sector. In order to better 
understand these, we construct the retarded propagator and discuss some of its 
properties in Section [S] It is then used to further analyse the planar divergences 
in Section |6l We conclude with a summary and an outlook. In Appendix El 
we discuss the relation to the matrix model introduced in [17] . Appendix [B] 
contains the calculation of the retarded propagator in momentum space. 

2 Setup 

We start by fixing some notation: For the commutation relations, we write 

[x'',x'']^ie'''' = iXl^e^"' (2.1) 

with 

■ 1 
-1 

Here we introduced a length scale Anc- The product of functions of these non- 
commuting coordinates can now be defined via the •-product, 

(/*.9)(^) = M^''''^5(^), (2.2) 

or by the twisted convolution 

(/ * gyCk) = {2n)-^ I d'k f{k)g{k - k)e'i'^^''"~'^" , (2.3) 

where the hat denotes Fourier transformation. For analytic functions, these are 
equivalent. When in doubt, we use (|2.3p . 



The Grosse-Wulkenhaar potential can now be introduced in the fohowing 
way: We define [HI IS] 

Here we defined 

A = f)-5A„c. (2.4) 

The choice il = 1, i.e., A = Anc, corresponds to the self-dual point. Obviously, 
and 



D^D^^' = -d^d^" T ^i\~'^e^^x''d^' - 4A -x„,r 



^-^.T^' 



The wave equation for a scalar field cf) of mass /i in the quadratic potential is 
then given by 

(-i (i?+i?+^ + D-D-^') + m') (^ = (a^a^ + A\-^x^x^' + I?) <p. (2.5) 

This is starting point for a discussion of the model in position space. The reader 
who is interested in this approach may thus directly jump to Section [5] The 
next two sections are devoted to the study of the model in terms of generalised 
eigenfunctions of the above wave operator for the case A — Anc- 



3 The eigenfunctions 

We want to study the eigenfunctions of the wave operator (12. 5p . For this, we 
restrict ourselves to the self-dual point, i.e., to the case A = Anc- The importance 
of this restriction is that then D^D^^ and Dj^D^'^ can be represented as •- 
multiplication from left, respectively right. In order to see this, we use the form 
(12. 2p of the Tk-product. For H — —-^x^x^, one obtains [17^, using e^rje = —rj, 

= -^ (x^x^ + i\^€,,,x^^d'' + iA^a^s^) / 



Analogously, one finds 



f^H^^D+D+f^f. 



2 -^M' 

Thus, if we find a complete set of orthonormal generalised eigenvectors \ks) of 
the Wigner transform of H, with eigenvalues k and a degeneracy index s, then 
we have 

(-1 {D+D+^ + D-D-^^) + n^) xti = (-A-^(fc + l) + /i^) Xl\. (3-1) 

where Xki i^ ^^"5 Weyl symbol of the ket-bra operator \ks){lt\. Furthermore, in 
this basis, the ^-product takes the form 

Xl\*xi^v^5{l~k')5''xl^v. (3.2) 



and because of the cyclicity of the integral we have 

d?x xi\ - 5{k - l)5^\ (3.3) 



Indeed, a basis with the required properties exists. As shown in [201 117) and 
below, the spectrum of H is the entire real line, with a two-fold degeneracy. 
The eigenvalues fc, I will be called the generalised momenta in the following. 

In order to find the eigenfunctions of H, we implement the commutation 
relations (j2.ip by choosing (we recall that here Anc = A) 

xo = Ag, xi = Ap, (3.4) 

where q and y = —idq are the position and momentum operators on L^(M). The 
Hamiltonian H thus becomes 

H = -^x^xl^ = 2(p2 -q^)^- {{q + p){q-p) + {q- p){q + p)) . 

We write this as 

H = -2(UV + VU), 

with 

U^^iq^p); V^^{q+p). (3.5) 

We have 

[[/, V] = i. 

Choosing the canonical representation for U and V, we thus have to solve the 
eigenvalue equation 

2i{udu + duu)ipkiu) = kipkiu), 

or 

ud^Mu) = (,-ij - ^) Mu) (3.6) 

Generalised eigenfunctions that solve this are given by [3T] 



1 

k 1 I 

2\/2tt "i 



V'^(w) = ^^U '^ 5 



"2 for w ^ 
otherwise 



It is straightforward to prove the orthonormality relations 

J2 [dkrkiu)rkiu') = 5{u~u') 

s '' 

There is a similar basis, obtained from ?/;^ by Fourier transformation, which is 
given by [H] 

4^N = geT^(-^*"^)r(-z| + i){u ± .6)'4-i 

The change back to the p, q representation is achieved by the unitary trans- 
formation 

Xiiq) = (2.^)-^ fdu e<--^^"+^)^±(u), (3.7a) 



riiiq) = (2.^)-^ fdu e<--^^"+-)4±(«). (3.7b) 



As shown in [20] , the results are the parabohc cyhnder functions that were used 
in [T7] and denoted by the same symbols (their convention is related to the one 
used here by £ = A~^fc/4). However, we note that U and V are multiples of the 
light cone coordinates, which are very convenient in two dimensions. Defining, 
cf. dMl) and ([33)1 . 



we obtain the Weyl symbol of the ket-bras |V'fc)(^*l in these coordinates as 

Xl\{u,v) = fdpe^P^{u-p/2\rk)mu + p/2). 



We compute this explicitly for the H — \- component. 

This vanishes for u <0. For u > 0, we obtain, using [521 (13-2.1)], 

xti'i^^y) = ^ / " dpe^p^iu^p/2)-^^-Hu+p/2y-^-^ 

^ "^ * ' dp e^'P'^^ {1/2 -p)-'^-i {1/2 +pyi^ 



47r j_i/2 
^ (2t.)'^ ^_,,„, r(»| + l)r(~»| + \)M{iL + 1, »1^ + 1, 4^H 

4^ "^ r(i^ + 1) 

where M is Kummer's confluent hypergeometric function of the first kind. For 

the component, one finds a similar expression, and also the H — and — h 

components can be expressed in terms of special functions, in this case Kum- 
mer's confluent hypergeometric function of the second kind. However, for the 
present discussion, the explicit form of the eigenfunctions in position space is 
not relevant. 



4 Quantisation in terms of the eigenfunctions 

We now want to discuss field theory at the self-dual point in terms of the con- 
tinuous set of generalised eigenfunctions. Expanding fields in terms of Xkn 



^Y. j '^^'^^ ^tixTk 



where 

we may write the wave operator W in matrix notation as 



with, cf. ([53]) . 



w^f /,;< = (-A-2(fc + /) + ^2) ^(fc _ k')s{i - ns^^'s'''. 



This can easily be inverted to yield a propagator 



ASttS 

^kl I'k' 



k + l-X^^i^ +ie(Jst{k,l) 



5{k - k')5{l - l')S"' 5* 



(4.1) 



Here astik,l) is a sign function which can be chosen such as to achieve the 
required causality properties for the propagator. Of course, the choice of this 
sign function affects the loop integrals we want to calculate later on. However, as 
we will see, some properties of the loop integrals are generic in that they do not 
depend on this choice. In particular, we find a divergence that is present even 
before evaluating the loop integral, and thus independent of the sign function. 
Hence, we will not invest too much care into a rigorous discussion of the possible 
choices for the sign function. 

For the graphical statement of naive Feyman rules, we use a double line 
notation, similar to [13j- The two lines can be interpreted as the bra and the 
ket of the eigenfunctions. According to (|4.ip . the propagator is given by 



It 



I't' 



Ktt' 



ks 



^ 6'' 6*U{k - k')Sil - I') 



-\^ 



k's' 



k+l-\^fi^+ie 



For the sake of notational simplicity, we dispensed with the sign function. How- 
ever, we have to keep in mind that the sign of e may depend on s, t, fc, I. 

From p.2p and p.3p it follows that in the scalar cj)^ model with coupling 
constant g, the vertex is given by 






l't_ 

ks 



= gS"" (5" (5"" S{k - k')S{l - l')S{j - f). 



k's' 



ju 



As a first application, we compute the fish graph in the (jr' model. One finds 
the following graphs: 




The first one is a planar and the second one a nonplanai|j graph. For the planar 

^In the notation of 1141 it has genus 3 = 0, but two boundary components (B = 2), and 
thus a hole. 



part, one obtains 

ig^X^S{k-k')S{l-l')S'''S*'' 



E/djd/ 



^ ^ [S{j~j')f- (4.2) 



k + j — A^/i^ + iei j' + I — X^ii'^ + i€2 



Here we introduced ei/2 in order to remind ourselves that the sign will in general 
depend on the generalised momenta and the degeneracy indices. Due to the 
presence of the square of a (5 distribution, this expression is divergent, already 
before evaluating the loop integral. It is thus no UV divergence in the usual 
sense (so also a cutoff in the spirit of [13] would not help). That is does not 
appear in the Euclidean version of the theory can be understood by noting that 
there, the loop integral is a sum of the form 

y^ S^f Sii' 1 1 . 

^r~f n + j + X^ii^ j' + 771 + \^p? 

Here, the square of the Kronecker 5 poses no problems. At the end of this and 
in the next two sections, we discuss the appearance of this divergence in some 
detail. But for the moment, we ignore it and formally absorb it in a divergent 
constant 5{Q) such that 5{i — j')^ = 5{Q)6{j — j'). We can then evaluate the 
loop integral in (j4.2p . and in particular discuss the unitarity of the model. We 
may write it as a convolution: 

^•?' ■ ^ . , ^ ■ ■ (4-3) 

J +ieik-l- J -162 

Using 

T[^J{p) = ^tV2^H{±p), (4.4) 

where J^ denotes the Fourier transform and H the Heaviside distribution, ancO 

Hi±p)H{±p) = Hi±p), 
H{±p)HiTp) = 0, 

we find that (|4.3p vanishes for ei = 62. In the case ei = —62, one obtains 

27ri- 



k — I + iei 



Because of ^^-pj^ ~ P^ ^ iTTS{x), the imaginary part of (j4.2p is thus given by a 
multiple of 

2Tr^g^X'^S{0)S{k - k')S{l ~ ly' S'*' S{k - I). (4.5) 

The multiplicity depends on how ei and £2 behave for the different combinations 
of s, t and u in (|4.2p . 

In the nonplanar graph all generalised momenta are fixed, so there is no loop 
integral to evaluate. We obtain 

txs't' 1 1 



k + V - A2/i2 + iej k' + 1 - A2,u2 _|_ j^^ ' 



ig'\'^6{k-l)5{k' -l')5''^5 

^Strictly speaking, these products are not well defined in the sense of Hormander's product 
of distributions |24| . Using Steinmann's concept of scaling degree |25l , one can show that the 
r.h.s. is, in a certain sense, the unambiguous extension of the product on the l.h.s. to the 
singularity at p = 0. 



Using again (|4.4I) . one can show that 



1 



1 



X zt ie X zt ie 



-dx 



1 



X ztie 



and that the products 



are not well-defined and have to be renor- 



x±ie x^ie 

maliseclJ. But even in that case, the product is well-defined on test functions 
vanishing in a neighborhood of the origin, so renormalisation ambiguities only 
affect the behavior at a; = 0. We may thus conclude that for k + k' ^ X'^^'^ the 
imaginary part of the nonplanar graph is given by 



g'X^S{k-l)S{k' -l')6''S 



st s:s' t' 



1 



(fc + fc' - A2/i2)2 ■ 



(4.6) 



For k + k' = A^^^ there are renormalisation ambiguities when ei = —62. But 
these are not relevant at the moment. The important point is that the contri- 
bution (|4.6p leads to a violation of unitarity. 

We now want to compute two graphs when the internal lines are put on 
the mass shell (and multiplied by 27r). For the planar graph, we again find 
the singularity due to the matching of generalised momenta at the two vertices. 
Writing this as (5(0) again, we obtain 



;«' 



2{2TTyg'X^SiO)Sik - k')5{l - l')6'' ,5" / dj 5{k +j- AV')'5(j + / - A^') 



8TT^g^X'^6{0)6{k - k')S{l - l')S''' S" 5{k - I). 



The factor 2 comes from the two-fold degeneracy. This is a multiple of (|4.5I) . 
Thus, it may be possible, by a suitable choice of the sign function, to fulfil the 
cutting rules. For the nonplanar part, however, we find 

47rVA^'5(fc - l)Sik' - l')S'''d''*' [d{k + k' - A^^)]^- 

We again find the renormalisation ambiguity at k + k' = X^^^, but no contri- 
bution of the form (j4.6p . Thus, unitarity is violated in a naive Feynman rules 
setting. 

Let us now come back to the subject of the strange divergence in the planar 
fish graph. It is straightforward to see that it is not specific to the cj)^ model, but 
also shows up in other planar graphs, such as the self-energy in the (j)'^ model, 
or the one-loop correction to the three-point function: 





^For a systematic treatment of renormalisation ambiguities in the products of distributions, 
we again refer to |25| . 



In fact, every ribbon graph in which a closed loop of a single line exists, 
i.e., every graph that contains a planar subgraph with a loop, is subject to this 
divergence. 

Furthermore, these divergences seem to be present also in a Yang-Feldman 
quantisation of the model: Then, one of the propagators in the loop integral is re- 
placed by a retarded propagator, and the other one by the Wightman two-point 
function of the free field [51 [TU] . The replacement of a Feynman propagator by 
a retarded one is unessential for the present discussion, since they should differ 
only in the ie description at the poles. Due to the broken translation invariance, 
there is some ambiguity in the definition of the free two-point function, but in 
any case it has to be a solution to the free field equation and it has to be com- 
patible with the commutator. If the retarded propagator can be written in the 
form (|4.ip . as we assumed above, then the commutator (which is derived from 
the retarded propagator) will conserve the generalised momenta. But then the 
two-point function must also have a component that conserves the generalised 
momenta. Thus, the strange planar divergences can not be avoided by using 
the Yang-Feldman formalism. In the following sections, we will show that they 
are no artefact of the use of an inappropriate basis, but also appear when the 
model is discussed in position space. 

Finally, we note that the planar divergences will also show up in the four- 
dimensional case. By a Lorentz transformation, one can always switch to a 
coordinate system where 6 is of the form 

Xfe 
XI, 

For the two spatial coordinates that now commute with time, the quadratic 
potential is the usual harmonic oscillator potential, so the generalised eigen- 
functions are given by 

where (f>mn are Weyl symbols of kct-bras of harmonic oscillator eigenstates. 
These fulfil 

Vklmn *Vk'l'm'n' = 0[l — k )0„m'0 Wkl'mn'^ 

instead of p.2p and p.3|) . The propagator will be of the form 

, . _i 

A 



klmn Vk'm'W ^_2(^ ^ ;) _ ^-2(^ + n) - ^,^ + Zea^Jfc, /) 

X 5{k - k')S{l - l')Srnn.'Snn'S'''S' 



■tt' 



Since the generalised momenta fc, I are conserved by the propagator and at the 
vertices, one will again find the square of a (5-distribution in the fish graph. 

In the setting of the generalised eigenfunctions, the planar divergences arise 
because the generalised momentum is conserved at the vertices and during prop- 
agation. One may thus suspect that the problem is absent when one is not at 
the self-dual point. In this situation, the description in terms of the generalised 



10 



eigenfunctions becomes quite complex. However, by switching to position space, 
one can show that the above reasoning is at least partially correct. As shown 
in the following sections, the singularity is absent (far enough) below the self- 
dual point, i.e., if the frequency of the potential is lower than the self-duality 
frequency, but present at and above the self-dual point. 

5 The retarded propagator 

In order to avoid the ambiguities in the definition of a Feynman propagatoip, we 
start by considering the retarded propagator. It can be constructed in position 
space, which avoids the use of the generalised eigenfunctions. In the massless 
case and without the quadratic potential, the wave operator for a scalar field is 
given by D = Adudy where we now use 

U ^ Xq ~ Xi, V = Xq + Xi, (5-1) 

instead of the notation p.Sp used in Section |31 The retarded propagator for this 
wave operator is 

Arot(wi,Wi;W2,W2) = \H(ui -U2)H{vi - W2), 

where H is again the Heaviside distribution. Thus, the square A^^^^ is well- 
defined without the need for any renormalisation (but see footnote [6]). In the 
presence of the quadratic potential, the retarded propagator will no longer be 
translation invariant, and the above propagator is multiplied with a function of 
ui,vitU2,V2. In the massless case, the wave operator for the Grosse-Wulkenhaar 
potential is given by, cf. (|2.5L 

4a„a„ -f- AX^^uv (5.2) 

and we have the following 

Proposition 5.1. The retarded propagator for the wave operator (J5.2I) is given 
by 



l^ret{ui,Vi\U2,V2) = \H{ui - U2)H [vi - V2)^{-lT ^^„,2» 



ul) {"l-vl) 



2"A2"n! 2"A2"n! 

n=0 



Proof. The series on the l.h.s. has infinite convergence radius and thus yields 
an analytic function T^(mi, wi;m2, W2)- This function (which is in fact a Bessel 
function) fulfils 

V{ui,vi]ui,vi) = 1, 
dmV{ui,vi\U2,vi) = 0, 
dv^V{ui,vi\ui,V2) = 0. 

The first equality assures that when both derivatives in duidu2 act on the Heavi- 
side distributions, then one still obtains a 6 distribution for coinciding points. 
Due to the other two equalities, the mixed terms, where one derivative acts on 



°It depends on the quantum state, which is not unique due to the lack of translation 
invariance. 
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a Heaviside distribution and the other one on V, vanish. Thus, it remains to 
show that 

{dmdvi + \^^UiVi)V{ui,Vi\U2,V2) =0, 

which is straightforward. D 

Before discussing the propagator in more detail, we express it in the coordi- 
nates 

Us^ui+U2, ut = ui- U2, (5.3a) 

Vs = Vi + V2, Vt^Vi- V2, (5.3b) 

as 

A,et(u„ .„.„«,) = ^2HMHMY.i-^r^^^^ (5.4) 

n=0 

= \H[ut)H{vt)Jo{\^'^ ^UtUsVtVs). (5.5) 

For imaginary arguments, i.e., for UsVg < 0, the Bessel function diverges as 
Jo{ix) ^ e^/-\/27ra; [231 (9.7.1)], which can be seen as the cause of the serious 
problems we will encounter. Using the inequality 

2\xy\<\x\^ + \y\^ (5.6) 

several times, one finds that the Bessel function is asymptotically bounded by 
a Gaussian, 

This is also true for the derivatives. It follows that the retarded propagator is 
well-defined on test functions that fulfil the bound 



with a = 47J37J7, where e can be chosen arbitrarily small. This is the Gelfand- 
Shilov space [25] Sc^a where a and A are the quadruples consisting of ^ and 
/^^ ; respectively, where e can be chosen arbitrarily small. For this, we will 
simply write 5q,^^(R"') with a — ^ and A = ^^ '^' in the following. By the 
above reasoning, Aiot can be interpreted as an element of S' a,A{^'^) 

Remark 5.2. We recall that for a massive theory (without quadratic potential) 
the retarded propagator is given by 

lH{ut)H{vt)Joi^J.^/utVt)■ 

One thus has the very natural interpretation of the propagator (j5.5p as the one 
for a position dependent mass //^ = X^'^UsVs- This is the value of the potential 
at the center of mass of the two points {ui,vi) and (1*2,^2). The problems for 
UsVs < stem from the fact that the model becomes tachyonic (and ever more 
so as UgVg — !> — cxo). 
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The Fourier transform of (j5.4p (which can be interpreted as an element of 
5'"' (M"*) with a and A as above) is 



,{KMMM) = -\Y.{^T b^-\K)8^-\ls) (rhiihr)^ ■ (5-7) 



Here fc^/t is the Fourier dual of Uj/t and l^jt that of v^jt- The appearance of 
derivatives of the (5-distribution in fc^ and Is indicates that momentum is not 
conserved. In Appendix |B] it is shown that the Fourier transform of (|5.5p can 
be expressed in terms of the Bessel function Kq. 

6 Planar divergences 

Having the retarded propagator at hand, we can now discuss the origin of the 
planar divergences found in SectionlH In the setting of the naive Feynman rules, 
the planar <^^ fish graph is given by 

where * denotes the •-product with Q replaced by —d. In the Yang-Feldman 
approach we would have to compute similar products, but with one of the 
propagators replaced by A^et and the other one by one of the two-point functions 
A-I-. This, however, requires the choice of a state, which we would like to avoid 
for as long as possible. We thus try to compute the product 

Arot(a;, y)*x*j,A,.ot(a;,2/). (6.1) 

Even though it has no direct physical significance, the study of this product 
helps to understand the origin of the planar divergences. We want to compute 
this product in the coordinates Us/t, «s/t- By using [m, u] = 2iA^(., we obtain 
the commutation relations 



{Us,Vs\ =[Mf,Wt] =[Ms,Mf] ^{Vs,Vt\ =0, {Us,Vt\ =["«<, Wg] = 4iA^ 



nc 



where [•, •] is the commutator where the *-commutator was used in the second 
argument. Thus, the correct twisting factor for our momenta is, cf. 



^-2i\l^(kJt+kth-'hkt-hk,} j-g 2) 

Now the question is in which sense the product (|6.ip should be defined. As al- 
ready noted, even the pointwise (commutative) product is not well-defined in the 
sense of Hormander's product of distributions. To our mind, the most conser- 
vative approach to a definition of (|6.ip is the followingj: In order to disentangle 
the problems connected to the ^-product and the distributional character of the 
retarded propagator, one begins by defining the •-product at different points. 
In the next step, one checks whether the limit of coinciding points makes sense. 
The definition of the ^-product at different points can be done by duality, as 
proposed in [27j : 

(A,et «)*,*„ Aret, I®g) = (Arot «) A,et , / ®*,5„ g) ■ (6.3) 



^In [9], this strategy was pursued for the definition of quasiplanar Wick products. 
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Here we wrote the planar •-product at different points in the form of a tensor 
product. Using (|6.2p . we have 






Formally, this may be written as 



f{ks,ls, kt,lt)g{ksjs,kt, h)- (6.4) 



(/®*,j„ g){us,Vs,ut,vt;Us,Vs,ut,vt) 



In order for the r.h.s. of (16.31) to be well-defined, we have to require the r.h.s. 
of (|6.4p to be an element of 5"^'^(R*) (or the r.h.s. of (j6.5p to be an element of 
5q.a(R*))- For this, we might have to choose / and g from a suitable subset of 
5Q^yi(K*). That this is possible if one is far enough below the self-dual point is 
the result of the following]J 

Proposition 6.1. For a = ^, A — /-'^ '^''^d, \fe\nc < (-^ — £). there is 
a nontrivial subset S of Sa,A{^'^), such that, for f,g ^ S, the r.h.s. of (16. 5p 
is well-defined as an element o/ 5q,.^(M^). More precisely, this is the case for 
S = 5f;I(M4) ^uh /3 ^ i, B = ^5^. 

Proof. According to [2S], the operator f{d) for an entire function / of order less 
than or equal to -g and type less than „i/g 2 is well-defined on the space 5^'^ 

and maps it to the space 5^ '^"^ . Using (J5.6I) . it is easy to see that the twisting 
in (j6.5p has order 2 and type A^^.. It follows that for /3 = ^ and B — —^ , 
the r.h.s. of (|6.5p is well-defined as an element of 5q,a'(IR^), provided that 
■\/eAnc < (A — e). It remains to show that the space S^'^{R^) is nontrivial. 
As can easily be seen by considering a Gaussian, this space, for a = /3 = i, is 
nontrivial provided that AB > 1/e, which is fulfilled. D 

By applying more sophisticated methods, it might be possible to get rid of 
the factor ^/e in the restriction on A. This would mean that to be below the 
self-dual point is a sufficient condition for the possibility to define the planar 
•-product at different points for elements of 5'q,/i(M^). But in any case we can 
show that it is a necessary condition: 

Proposition 6.2. For A„c > A, a = ^ and A = r^ ; there are no e,e' > 
such that there are nontrivial f,g ^ S°'' (M*), with A' = ^ ^'^ ' for which 
the r.h.s. of ^^ is an element ofS^^^iR^). 

Proof. We assume that such e,e' and such functions /, g exist. Now, according 
to [25], elements of 5^'^, for /3 < 1, are entire functions that fulfil the bound 



-'Similar considerations on the ^-product of elements of 5 '^ can be found in I28| . 
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with b = — -{Be)^^^ + S, where S can be chosen arbitrarily small. In our case, 
this means that / and g are entire functions that fulfil the bounds 

\fix + iy)\<ce'''\y\\ \g{x + iy)\<c'e'''\y\\ (6.6) 

with 6' = (A — e')^ + 6' . Here x and y are elements of M*. Since we assumed 
that the r.h.s. of (16. 4p is an element of 5"'"^(R^), also the bound 



|(/ ®.^,„ gYix + ty-S: + ty)\ < CeKl^l'+IS^, 

with b — {X — e)-^ + 5, has to be fulfilled. We define the matrix 

/O -1\ 
10 
0-100 

yi y 

The above inequality then leads to 

e'^U\-\'-\yn\f{z)g{^Jz)\<Ce''i^y\'+^-n, 

where z = x + iy. Thus, we have 

\f{z)gitjz)\ < Ce-(2^"=-'')l-l'+(2^-+'')l^l'. (6.7) 

Because of Anc > A, (5 can be chosen such that 2Xf-^^ — 6 is positive. Thus, the 
function f{z)g{ijz) falls off with order 2 and type 2X^^ — b in the real direction. 
On the other hand, from the bounds (J6.6I) on / and g, it follows that 

\f{z)g{ijz)\ <C'e'''(l"l'+l'yl'). (6.8) 

Thus, f{z)g(i^z) has growth of order 2 and type b' in the imaginary direction. 
Hence, for 2X^^ — b > b' , the entire function F{z) = e^ ^ f{z)g{ijz) is bounded 
on the real and imaginary axes and, by (j6.8p . uniformly bounded by C'e^^ '^' in 
between the axes. By the Phragmen-Lindelof principle [IHl Thm 2.5.2], it is thus 
bounded on the whole complex plane, and can only be a constant. But since its 
limit in the real direction is 0, cf. (|6.7p . we have F{z) = 0. Thus, nontrivial /, g 
canonlyexistfor 2A^c-fe< b' . This means 2A^^-(A-£)2 -(A-e')^ < 6 + 6'. It 
is clear that for Anc > A, ^ and 6' can not be chosen arbitrarily small, contrary 
to our assumption. D 

By this analysis, it is not possible to define the planar product of two re- 
tarded propagators at different points at or above the self-dual point. Thus, the 
problem in defining the product (|6.ip does not stem from the limit of coinciding 
points, since already the product at different points is ill-defined. In this sense, 
this is no UV divergence. 

The result of this rather abstract argument can be checked with a concrete 
calculation. We will do this in a formal way, i.e., we use the series form (|5.4p (or 
equivalently (15. 7p ) of the retarded propagator and compute the planar product 
for the individual terms. The hope is that we obtain a (power) series that can 
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again be summed up. The twisted convolution that we want to compute is 
theiS 






mn 

m+1 



^ X! / dfcsd/sdfctd/ 



n+1 



We now use 



<5(")(x-y)/(a;) = f](-ir(^)/(™)(y)<5("-'")(:E-y). 

m=0 ^ ^ 

to get rid of the twisting factor, by first applying this equality to kg and Is and 
then to ks and Is- One thus obtains 



, "-J4 + 1 / \ n-J3 + l 



kt — kt—iej \lt—lt—if- 

Let us consider the integral over kt. For m — J2 > and n — J4 > 0, the integral 
yields, in position space, a multiple of 

For m — J2 < one has n — J4>J2~m>0 and the integral yields, in position 
space, a multiple of the product of 6^^^~"^~^\ut) and H{ut)ul~^* ■ Albeit this 
product is not well-defined in the sense of Hormander, it vanishes in the sense 
of Steinmann's scaling degree. The analogous argument works for n — J4 < 0. 
We thus obtain, in position space, a multiple of 

00 min(rn,n) 3 -^J2j^ 

^nlii \'R(i^\ \^ \^ n\n\m\m\ l^^nc^ / 1 yi+J4 

m,n^O ji—0 

Apart from the factor (— ipi+^i the summand is invariant under the exchanges 
ji '^ 33 a-iid J2 <->■ ji- It follows that only terms where J2 + J4 and ji + J3 are 
even contribute. We can thus write the above as 

00 min(m,n) 4 j._^; 

m,n=0 kd=0 ^ ' 



(6.9) 



^^ Since our calculation is formal anyway, we could also compute the ^-product with the 
formal series 112.211 . The result is the same. But since a calculation in momentum space is 
needed later on, we chose to present it in terms of the twisted convolution 112.31 1. 

16 



with 

min(2fc,min(m,n)) 



mn 



n\m\ /'_-| \j 

i!(2fe-j)!(m-j)!(m-2fc+j)!(n-j)!(n-2fe+j)! V ^1 ■ 
J— max(0,2fc— min(m,n)) 



For m — n — k, one finds c™^ = (—1)'". While it seems to be hard make a 
statement on the convergence of the series in (I6.9P for fixed but general utVg 
and UgVt, it is easy to show that it does not converge for Ug — Vg — 0. This 
also shows that it does not converge as a power series, contrary to (|5.4L since 
the zeroth order coefficient does not converge. For Us = u^ = we only get a 
contribution for ttt, = n = A; = ^, so that the above series reduces td^n 

E(^) • (6-10) 

m=0 ^ ^ 

Obvioulsy, this diverges unless Anc < A, so the planar product (|6.ip is only 
well-defined below the self-dual point. Using (|2.4p . the above reduces to 

oo _, 

X ^ ^/l^ t 



n" 



^^ 1-^4 

m=0 

For i7 = 1 — e, wc thus find a divergence e^^ as e — > 0. Note that the problem is 
not that the loop integral over the momenta diverges. Thus, this is no ordinary 
UV divergence, similarly to what we found previously. 

Remark 6.3. Let us consider what happens in the case of a massive field. Then 
the series in (|5.4I) will be a power series in A^^" and /i^. However, at zeroth 
order in ^^, one finds again (J5.4I) . Thus, when one calculates the planar product 
in a formal way, then a nonvanishing mass does not help. If one does not resort 
to a formal calculation, then it is to be expected that the propagator is still only 
defined as an element of 5'q.a(K.4) with a and A as above, since in the long 
range the quadratic potential will always dominate the mass, cf. Remark 15.21 

The argument given up to now is not complete in the sense that we computed 
the planar square (j6.ll) of the retarded propagator, which is not what appears in 
actual loop calculations. In the Yang-Feldman formalism on noncommutative 
spacetimes, the planar fish graph loop integral is of the form \Q_, J^ 

A+(.T,y)*:j*yAret(2;, y) + A_(a;,y)*:r*y \ct{x,y), (6.11) 

where A_|_ is the Wightman two-point function and A_(a;, y) — A+(y, x). Thus, 
one has to choose a state, which, however, is not unique due to the lack of 
translation invariance. But usually the two-point function is defined on the 
same test function space as the retarded propagator (or a subset thereof) . Thus, 
by the analysis in the beginning of this section (which only used the structure 
of the test function space), we expect the same problems as above. In order 
to be more concrete, we choose a particular two-point function and repeat the 
(formal) calculation from above. 
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For the coefBcient of the utUsVtVs component, one finds, using c™^ , j = {m + 1)(— 1) 
the series 5I]m("^ + l)^(-^nc/A)*"', which diverges even worse. 
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A two-point function A-|_ has to be compatible with the commutator, which 
is defined via the retarded propagator, i.e.. 



A+{ks,ls,kt,lt) - A+{ks,ls,-kt,-lt) 

= i [A,.ct{ks,ls, ktJt) ~ Ai.ct{ks,ls,-kt,-k) 



Furthermore, it has to be a solution of the field equation. Finally, some kind of 
positivity would be nice. Thus, ignoring the usual infrared problema^'^l. 



n=0 






would be a suitable two-point function, cf. (|5.7p . Using this two-point function, 
we compute the product (16.111) at Us = Vs = 0, i.e., we consider the component 
where all derivatives of the S distributions of kg, Is, kg and Ig are shifted on the 
twisting. For given m > n we obtain 



(47r)^(m-n)! \ X ) 



ihY, 



dkfdlt 



n— 771+1 






Here we used 



c"5(")(a;) = (-l)"n!(™j5("'~")(x). 



Shifting the derivatives w.r.t. It and fcf away from the 5 distributions, we see that 
this vanishes for m > n. For m < n, we also get a vanishing expression, since it 
involves the products l^5'^"^\lt) and k^5'^™\kt). Thus, only the contributions 
with m = n survive. As above, these are independent of rm^ so we again find 
the series (|6.10p . which diverges for Anc > A. 

Remark 6.4. Even if one is not at the self-dual point, one still has, under a suit- 
able exchange of positions and momenta, the duality S[4>, m, Q] — >■ Q,'^S[(I>, mi}^^, ft' 
cf. [TS]. Thus, one might wonder about the compatibility of this fact with the 
above finding that the model behaves well for ft < 1 but diverges otherwise. 
The point is that the above duality does not respect crucial properties of quan- 
tum field theories on Minkowski space, such as causality or positive energy. In 
particular, the retarded propagator is not invariant under the above duality 
transformation. 

Remark 6.5. Also on Euclidean space divergences at the self-dual point were 
found, namely in planar tadpoles of the Gross- Neveu model [301 App. A. 4]. 
Similar to the findings presented here, the singularity is present even before 
the loop integral is evaluated. Also there, the origin is the behaviour of the 

^^By restricting to test functions that vanish in a neighborhood of kt = It = 0. 
^^The integral itself is UV finite, but has the usual infrared problems. 
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propagator for large spatial distancea^H It would be interesting to further study 
the similarities of the two effects. 

Remark 6.6. Finally, a comment on the (j)'^ model. The fish graph loop calculated 
here also occurs in the four-point function of the (j)^ model, so the problem is not 
specific to the (p^ model. For the two- loop self-energy graph shown in Section 3] 
one has to compute products involving a retarded propagator and two of the 
two-point functions A±. If already the product of Aret with one of these does 
not exist, then neither do the higher order products. 

7 Summary &; Outlook 

We discussed noncommutative field theory with Grosse-Wulkenhaar potential on 
the the two-dimensional Minkowski space in two ways: In the first approach, we 
restricted ourselves to the self-dual point and used a continuous set of generalised 
eigenfunctions of the wave operator. This we used to postulate naive Feynman 
rules. In this setting, we found a new type of divergence in the planar sector. 
By considering the situation in position space, we showed that this divergence 
is not due to an inappropriate choice of the basis. Instead, the fast growth of 
the propagator in some directions makes the definition of the planar -A-product 
impossible, even before considering the limit of coinciding points. 

In our opinion, the appearance of this new type of divergences is an inter- 
esting phenomenon that deserves more detailed studies. These could proceed 
along the following lines: In order to relate the two approaches discussed here 
(position space and generalised momenta), it would be useful to have a repre- 
sentation of the retarded propagator in terms of the eigenfunctions Xm- This 
would amount to find an appropriate sign function ast{k,l)- 

Furthermore, it would be interesting to know whether one can get rid of the 
factor ^/e in the restriction on A in Proposition 16. II This would mean that the 
model is well-defined on the whole interval Vt e [0, 1). Otherwise, the self-dual 
point may not be so special after all. Another (possibly related) question is the 
following: We have shown that at and above the self-dual point the individual 
terms of a series expansion of the planar square of the retarded propagator 
diverge. We conjecture that below the self-dual point all indivual terms in this 
expansion converge. If this is the case, it remains to check whether the series as 
a whole converges below the self-dual point. 

One could also study the model in the matrix basis of Grosse and Wulken- 
haar. In the Minkowski case, the propagator will then take a more complicated 
form than in the Euclidean case, but it might still be possible check whether 
one runs into problems similar to those discussed here. 

Another important point is renormalisation. Because of the uncommon type 
of the divergences, it is not clear whether such a program can be successful 
and how one should proceed, but perhaps the formal renormalisation used in 
Section 13] would be a good starting point. 

Finally, one should treat the four-dimensional case. As discussed in Section 
m the use of the generalised eigenfunctions will again lead to divergences in 
planar (sub)graphs. Preliminary results suggest that this is also true in posi- 



^^The author would like to thank F. Vignes-Tourneret for private communication on his 
work. He would also like to thank the referee for pointing out this reference. 
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tion space, at least in a formal sense. In that case, it would be important to 
understand why these problems are absent in the Euclidean setting. 
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A The relation to the matrix model 

We now want to clarify the connection to the matrix model setting proposed by 
Fischer and Szabo [T7]. They work at the self-dual point and also consider the 
eigenfunctions x±, ry±, cf. (|3.7|) . However, they suggest to transform the model 
to matrix form by considering the Gelfand triple 

5^(M)cL2(M)c5^(K)', 

where 5"(R) is a Gelfand-Shilov space [55]. For elements of this space, they 
claim (their Theorem 4.2), 

hm {xi\4>)=Q yc^eS^iR), 

k-^—ioo 

, lim {v'±m^O V0e5^(i?), 

where the limit is taken in the lower (upper) complex half-plane. Furthermore, 
the eigenfunctions x± and 77^ have poles at fc = —2i{2n+l) and k = +2i{2n+l), 
respectively, for n £ Nq. The corresponding residues are given by 



Res [xi(<z)] « f-iq) cc e+'-Hr^iV^iq), 

k=-2i(2n+l) 



Res Ml)] « fnil) « e '-HnW+iq), 

k—-\-2i{2n+l) 

where Hn are the Hermite polynomials. Finally, one has 

xtiq) ^Vkii)- 

From these facts they conclude (their Corollary 4.3), that, by closing the contour 
of integration^^ 



= J2jdk\x':){x'M = J2\fn){fn 



Since the /~ are neither elements of 5"(]R) nor of L'^{R), the convergence seems 
to be in iS"(]R)', but it is not clear in which topology. 

Now the (f) in the above equation is still only a ket. However, for ip = Icf)) {ipl e 
iS2(M) (g) 5^(R), one obtains the expansion 



-"^^To be precise, in the mentioned Corollary, they write (j) = o Es Enl/" '')(/nl</')- The 
second term seems to be added for convenience. 
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with 

¥'m„=(/+l0)(V^I/n)eC 

and 

Jmn — \jm)\Jn \- 

The fmn fulfil the usual properties of a matrix base, i.e.r^l 

jTnnJiyi'n' — '^nm' Jmn' : ^^ Jmn — ^mn- 

One can thus use them to bring the model into matrix form and treat it similarly 
to [13]. One then arrives at the following representation of the propagatoq^ 
(equation (3.58) of 17 with a — i): 

^mn n'm' — ^mm'^nn' TTl ; ' TT ; 7~^ ^ l-^'-^j 

—ii[m + n + 1) + X^ii^ 

However, some remarks are in order: As mentioned in [17, the use of 
Gelfand-Shilov spaces as test function spaces for noncommutative field theo- 
ries has been proposed by several authors [2H1 EZ] • This would imply that fields 
are elements of the dual space 5^(M'^)'. In the setting of [17], however, the 
fields are elements of the Gelfand-Shilov spaces iS"(M^). Thus, the fields are 
vanishing rapidly at infinity. If the quadratic potential was absent, it would be 
clear that this would not be a suitable space for the fields to live in, since it 
would contain no solution of the free equation of motion. But the quadratic 
potential does not change this, as can be seen from the absence of poles in the 
propagator (jA.ip . Thus, the space of fields proposed by [17] does not contain 
the solutions of the free field equation, even though such solutions exists, as 
can be seen from the pole in (|4.ip . To disregard the solutions of the free field 
equation is certainly a deviation from the principles of perturbative QFT. In 
particular, it is not clear how to describe asymptotic states (and thus to allow 
for a contact with experiment). 

But also from a mathematical point of view, the approach followed in |17| 
seems to be questionable, as their basic Theorem 4.2 is incorrect. This can be 
seen by the following counterexample. In Theorem 4.2 it is stated, that 

lim {xl\4>)=(i V0e5:(R). 

fc— f — 'ioo 

In order to test this assertion, we choose 0(q) = e^"'^ "^ , with some real 
constant a. We have XkilY — Vkil) ^^"^ (^^^ conventions used here are related 
to those used in [17 by f = A^^/c/4, £" ^ A^^, i^ — — x~5' ^^^^ *^^ parameters 
of dl] on the l.h.s.) 



V^{q) = az-Vr(f + ^)D_^.{TV2lX-^q), 



^^However, it is not clear in which sense these relations should be understood, since the 
/„ are not in L-^(R). Also an interpretation in a distributional sense as for I I3.2I I. I I3.3I I is not 
possible, since one can not interpret /„ as a distribution in n, due to the discreteness of the 
imaginary eigenvalues. 

-"^^Using the basis r; as a starting point, one would arrive at a similar propagator where the 
denominator is replaced by its complex conjugate. 
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where C is some constant independent of k and Di, is a parabolic cylinder 
function. For ^tv < 0, it is given by 



1 /"Oo 

i 1 ,2 



For "^v < 0, we compute 

/oo 
-CX3 

^ ^14 / aq e-^^e-V / di e±^'*e-Vi- 

r(-i^) 7-00 Jo 

= — / di / dg e —^i'^2^l e^^e~^t 

r(-t^) Jo J-oo 



A2-5 / 47r \5 /-"^ , t^,, 2> ^ ,^ , 

d-U e 2 2a + .u 2 



Ar(-|) / TT V /12a-i 



Here we supposed that ^f^f^qif > 0, i.e., a > ^. It follows that 
It is obvious that this does not converge for k — ^ — ioo. 

B The retarded propagator in momentum space 

We compute the Fourier transform 

Arotl^sj is,kt,lt) 

of the retarded propagator. From Remark 1 5. 2 1 we know that the retarded propa- 
gator can be interpreted as the one for a position dependent mass. The Fourier 
transform w.r.t. kt and It is thus well-known, and we obtain 



We now consider the cases UsVg > and UgVg < seperately. In the first case, 
we use the coordinates x — y/u/ul, y — y^Ug/vs and obtain the integral 



|.oo POO 2x 

j^ / Ax dy — cos{ksXy + ^-f) 
•/ i' y 



« {kt-ie){lt-ie)~f^, 
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In the case UsVs < we instead find 

vrrrj / da; / dw — cosiksxy — —) 5- 

Now for a, 6 > one lias [31, (3.868.2 & 4)] 







dx cos(a^x+ ■^)i = -7ryo(2a6), 



/•OO 

/ dx cos(a2a; - ^)i = 2A'o(2a6). 
Jo 



/o 
In the case sign ks — sign Is we thus obtain 



-4- / AxY^{2x^\kj:\)- — ^ 

/o (fct -2e)(/t-ze)- ^ 



(2t) 



8 



dx Ko{2xy/\kJ7\)- 



^ ^ Jo (fct-ie)(/t -ie) + ^ 

and in the case sign ks 7^ sign Is we get 



^ / dx ii:o(2a;v/|fcX|) 



{kt - ie){lt - i€) - ^ 



+ j^ dxYoi2x,/\kj7\)- — --^. 

^ ' Jo [kt - ie){lt - le) + jj^ 

We have the asymptotic relation [23, (9.7.2)] 

and may thus change the contour of integration for the integrals involving Kq 
to — > — loo. In the first case, we pick up a pole li kt + k < and in the second 
one if kt + lt> 0. We may then use [H (9.6.5)] 

nYoiiz) = nrloiz) - 2Koiz) 

and Iq{z) = Jo{iz) to obtain 



^AX^Hi-kt ~ lt)H{ktlt)Ko{-4iX-'^\kslsktk\ + e) 

S-rri A\4: 



+ ^AX'^Hi^h - lt)H{-ktlt)Koi'iX'y/\kslsktlt\- le) 



-Atzi 



00 



dx Jo(2xV|^X|)- 



^ ' Jo f^ - {kt - ie){lt ~ le) 

for sign ks — sign Ig and 

^4X''H{kt + lt)H{ktlt)Ko{^X'^\kslsktlt\ - te) 

+ 0^AX''H{kt + lt)H{-ktlt)Ko{~^iX^V\kslsktlt\ + e) 



(2^)" 



X 



dx Jo{2xy/\ksls\)- 







■§j + {kt -ie)(lt -ie) 
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for signks ^ signZs. We have [ST, (6.532.4)] 

X 

dx — —Jf)(ax) — KQ{ak), a > 0, 5Rfe > 0. 



Q X -r ft^ 



In order to solve the above integrals, we thus have to choose the root with 
positive real part of ^{kt — ie){lt — Je). For signfcg = signZs, we obtain 



Stti 



(2,r) 



■^4X''H{~kt - lt)H{ktlt)Koi~4iX-'^\kJsktlt\ + e) 



+ ^AX^Hi-h - lt)H{-ktlt)KoiAX'y/\kJsktlt\- ze) 
+ ^AX''H{ktlt)Ko{sign{kt + lt)4iX'^\kJMt\ + e) 

+ 0^AX''H{-ktlt)K^{AX^^%lMt\+iesign{kt+lt)). 
This can be written as 

^AX" sigu{kt + k) (^H{ktlt)Ko{sign{kt + lt)AiX'^\kJMt\ + e) 



+H{-ktlt)Ko{4:X'^ y/\kJMt\ + iesigii{kt + k)) 
For signfcs ^ sign/^, one likewise obtains 

^AXUignih + It) (^H{ktlt)KoiAXW\kJsktk\-iesign{kt+lt)) 

+H{-ktlt)Ko{~ sign{kt + ItHiX^y/lkslshkl + e)) . 
In total, we thus obtain 

ArctiksJs, kt,k) = -■(l^4A'* sign(A:t + k) iiign{ksls) 

X (^H{ktltksls)Ko{s\gn{kt + lt)s\gn{ksls)^iX^ \/\kJM^\ + e) 



-H{-ktltksls)Kn{AX^ y^\kj,ktlt\ + ie sign(fct + k) sign{ksh)) . 



Note that for large kgu,lsu, this is bounded, but highly oscillatory in some 
directions, as expected for an element of S'"' (W^) with a — ^. 
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